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Abstract. Discrete bright breathers are well known phenomena. They are 
localized excitations that consist of a few excited oscillators in a lattice and 
the rest of them having very small amplitude or none. In this paper we are 
interested in the opposite kind of localization, or discrete dark breathers, where 
most of the oscillators are excited and one or a few units of them have very 
small amplitude. We investigate using band analysis, Klein-Gordon lattices at 
frequencies not close to the linear ones. Dark breathers at low coupling are 
shown to be stable for Klein— Gordon chains with soft on-site potentials and 
repulsive dispersive interaction, and with hard on-site potentials and attractive 
dispersive interactions. At higher coupling dark breathers lose their stability 
via subharmonic, harmonic or oscillatory bifurcations, depending on the model. 
However, most of these bifurcations are harmless in the sense that they preserve 
dark localization. None of these bifurcations disappear when the system is infinite. 
Dark breathers in dissipative systems are found to be stable for both kinds of 
dispersive interaction. 
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1. Introduction 

It is well known that intrinsic localized modes (also called discrete breathers) are exact, 
periodic and localized solutions that can be obtained in a large variety of nonlinear 
discrete systems. They are becoming a new paradigm for understanding many aspects 
of the behaviour of discrete systems (for a review, see, e.g., ||). Mackay and Aubry 
1^ proved analytically their existence and the conditions for their stability [Q , under 
rather general hypotheses. Since then, many accurate numerical methods have been 
used to obtain breathers as exact numerical solutions up to machine precision 
which permits the analysis of breathers properties. Thus, for a given model, it is 
possible to perform a numerical study of the ranges of existence and stability in the 
parameters space. 

The term discrete breather is usually understood as a localized, periodic solution 
in a discrete system, but with a small number of excited oscillators. When only one 
oscillator has large amplitude it is known as one-site breather. When more than one 
oscillator have large amplitude, the term multibreather is used. ^Hereafter, we will 
use the term bright breather when there are one or a few oscillators vibrating with 
large amplitude whereas the rest of them oscillate with small amplitude. 

However, localization can be manifested in a different way, which consists of all 
the oscillators vibrating with large amplitude except one or a few of them oscillating 
with very small amplitude. The natural name for these entities is dark breathers., 
analogously to the well known term dark solitons. For the nonlinear Schrodinger 
equations, which governs both nonlinear optical modes in fibers and dilute Bose- 
Einstein condensates, two different kinds of scalar soliton solutions, bright and dark, 
are known ^. Thus, a dark soliton is a solution which has a point with zero 
amplitude, that is, a soliton defining the absence of matter or energy. Since then, 
many papers have appeared referring to theoretical and experimental results referred 
to this entities. The effects of discreteness on the properties and propagation dynamics 
of dark solitons has been analyzed in the discrete nonlinear Schrodinger equation 
(DNLS), which is thought to be a good approximation for frequencies close to the 
linear frequency 0, ^, y, the last one, also with numerics on actual Klein- 

Gordon systems. It is worth remarking that some examples of dark localization have 
been observed experimentally , and some structural properties have been analyzed 
in(l|. 

In this paper, we perform an analysis of the existence and stability of dark 
breathers for different models based on the properties of the band structure of the 
Newton operator. We have found that there exist stable dark breathers for a variety 
of one-dimensional Klein-Gordon lattices. For soft on-site potentials, dark breathers 
are stable only for repulsive dispersive interaction and for hard on-site potentials 
the stability founded is for attractive interaction. This results agree with the ones 
that have been found for the DNLS approximation in as a consequence of the 
modulational instability of the constant amplitude background. 

It is not clear in which physical systems dark breathers can play a significant role. 
A possibility is DNA, we conjecture that dark breathers can occur in DNA chains at 
high temperature, close to thermal denaturation JTsf . In this situation, a great number 
of molecules are vibrating with high amplitude, whereas a few of them could be almost 
at rest. 

This paper is organized as follows: In section ^, we describe the proposed Klein- 
Gordon models and the resulting evolution equations. In Section |^, we expose the tools 
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for calculating dark breathers, and we explicitly show that the theorem of existence 
of Mackay and Aubry gives an afRrmative answer to the question of existence of dark 
breathers. In section ^ we expose the method for the analysis of the linear stability 
of breather solutions using both the Floquet multipliers and Aubry band's theory. 
In Section ^, we investigate the stability of dark breathers for chains with soft on- 
site potential and attractive interactions between the particles. We show that it is not 
possible to obtain stable dark breathers for every value of the coupling parameter. This 
negative result suggested us to consider the study of chains with repulsive interactions 
between particles. Our results confirm that in this case there exist stable dark 
breathers up to significant values of the coupling. In Section ^, lattices with hard on- 
site potentials are considered. For these cases, dark breathers are stable provided that 
particles interact through an attractive potential. Dark breathers become unstable 
through subharmonic, harmonic and oscillatory bifurcations depending on the type of 
the on-site potential. In section |^, we show that both subharmonic and oscillatory 
bifurcations do not disappear in infinite systems. The paper concludes with a short 
summary of the main results and some prospective in Section ^. 



2. Models 



We study one-dimensional, anharmonic, Hamiltonian lattices of the Klein-Gordon 
type. The Hamiltonian is given by 

H = J2ilul + VM) + eW{u) (1) 

n 

where m„ are the coordinates of the oscillators referred to their equilibrium positions; 
V{un) represents the on-site potential; u represents the set of variables {u„}; and 
eW{u) represents the couphng potential, with e being a parameter that describes the 
strength of the coupling. We suppose initially that e is positive and W{u) is given by 

n 

This interaction is attractive because a nonzero value of a variable tends to increase 
the values of the neighbouring variables with the same sign. The on-site potential is 
given by 

V{Un)^^Ujlul+ (j){u„) (3) 

with (t>(u„) being the anharmonic part of the potential. The variables are scaled so 
that all the particles in the lattice have mass unity and the linear frequency ujo = 1. 
The dynamical equations for this system are 

iln + uj^Un + cj)'{un) + e{2u„ - w„_i - 1i„+l) = 0. (4) 

These equations do not have analytical solutions and must be solved numerically. 
The solutions depend obviously on the chosen potentials V{un) and W{u). We will 
analyze the system for several V{un) and coupling interactions W{u). These models 
appear in many physical systems, a known example, with a suitable on-site potential 
V{un) is the Peyrard-Bishop model for DNA where the variables m„ represent 
the stretching of the base pairs. 



Dark breathers stability 



4 



3. Dark breathers existence 

We look for spatially localized, time-reversible and time-periodic solutions of 
equations (||) with a given frequency uit and a continuous second derivative. Therefore, 
the functions u„(i) can be obtained up to machine precision by truncated Fourier series 
of the form 

Unit) = zo+ 2z^cos(fcwbt). (5) 

k=l 

We distinguish three types of solutions of the isolated oscillators that can be 
coded in the following way: an — for an oscillator at rest {un{t) = 0, Vt); (t„ ~ +1 
for an excited oscillator with frequency Wb and u„(0) > 0; finally, cr„ = — 1 for an 
excited oscillator with frequency Wb and w„(0) < 0. Time-reversible solutions of 
the whole system at e = (anticontinuous limit), can be referred to by a coding 
sequence a = {<Tn}- Therefore, a = {0,...,0,1,0,...,0} corresponds to a one-site 
breather. Other codes can be cr = {0, 0, 1, 1, 0, 0} for a symmetric two-site 
breather, and a — {1, 1, 0, 1, 1} for a one-site dark breather (this case correspond 
to the background in phase). 

The method for calculating dark breathers is based on the general methods for 
obtaining breathers g |l^, |l§. 

The existence theorem by MacKay and Aubry Q establishes that every solution 
at the anticontinuous limit, corresponding to a code sequence can be continued up to a 
certain value of the coupling parameter 7^ 0, as long as the following two hypotheses 
are fulfilled: 

• The orbits of the uncoupled excited oscillators with the chosen frequency have to 
be such that ^ 0, where I = J pdq is the action variable of the oscillator. 
That is, the oscillator is truly nonlinear at that frequency. 

• The frequency of the orbit must be such that pcui, ^ loq for any integer p. That 
is, none of the breather harmonics coincide with the linear frequency loq. 

Therefore, this theorem gives an immediate answer to the question of the existence 
of dark breathers as they are obtained by continuation of the configuration mentioned 
above a = {1,...,1,0,1,...,1} . However, the theorem does not give an estimate for 
the value of the coupling Sc where the dark breather ceases to exist. Dark breathers 
have to be calculated numerically for each value of the coupling parameter, also it is 
worth investigating whether dark breathers are stable or not, and if they are stable, 
up to which value of the coupling parameter. 

4. Dark breathers stability 

The stability analysis of a given breather solution can be performed numerically 
[0, |2^, |2l) . The linearized equations corresponding to perturbations of this solution 
are 

^r, + Wo in + (t>"iUn) ^„ + £(2^„ - ^n-l - Cn+l) = (6) 

where ^ — {f„(i)} represents a small perturbation of the solution of the dynamical 
equations, u{t) = {u„(t)}. The linear stability of these solutions can be studied 
by finding the eigenvalues of the Floquet matrix To, called Floquet multipliers. 
The Floquet matrix transforms the column matrix with elements given by Cn(0) 
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and 7r(0) = ^n(O) into the corresponding column matrix with elements £,niTb) and 
'n{Ti,) = ^niTi,) for n — 1... m and T(, — 2tt/uji,, that is 



{7r„(T,)} J-'"'\ K(0)} 

The Floquet matrix J-q can be obtained numerically by integrating equations a 
breather period. In order to get accurate results [ p2[ , we have used a symplectic 
integrator. Equation (^) can be written as an eigenvalue equation 

iN-iu{t),e)-On ^ E^n (8) 
where Af{u{t), e)} is called the Newton operator. The solutions of equation (H) can be 
described as the eigenfunctions of M for E — 0. The fact that the linearized system 
is Hamiltonian and real implies that the Floquet operator is a real and symplectic 
operator. The consequence is that if A is an eigenvalue, then 1/A, A* and 1/A* are 
also eigenvalues, and therefore a necessary condition for linear stability is that every 
eigenvalue has modulus one, that is, they are located at the unit circle in the complex 
plane. Besides, there is always a double eigenvalue at 1 + Oi from the fact that the 
derivative {un{t)} is always a solution of (||) with E = Q. 

The Floquet multipliers at the anticontinuous limit (e — 0) can be easily obtained 
for bright and dark breathers. If an oscillator at rest is considered, equation 
becomes 

L+^Un=0 (9) 

with solution (^{t) = (.0^^'^°*, and therefore, the corresponding eigenvalue of is 
A = exp(i27ra;o/w6)- If an oscillator is excited, equation (||) becomes 

^« + Wo ^« +(/)"(«„) ^„ =0. (10) 

This equation has Un{t) as a solution, which is periodic and therefore with Floquet 
multiplier A = 1. Thus, for a bright breather, we have, taking into account 
their multiplicity, 2(A^ — 1) eigenvalues corresponding to the rest oscillators at 
exp(±i27rtJo/'-^&) and two at +1 corresponding to the excited one. For a dark breather 
at the anticontinuous limit there are 2{N — 1) eigenvalues at +1, and a couple of 
conjugate eigenvalues at exp(±i27rijJo/a;b). 

Figure |l| shows the Floquet multipliers for both a bright and a dark breather at 
the anti-continuous limit e = 0. When the coupling e is switched on, these eigenvalues 
move on the complex plane as continuous functions of e, and an instability can be 
produced only in three different ways : (a) a couple of conjugate eigenvalues reaches 
the value 1 + Oi {9 — 0) and leaves the unit circle along the real line (harmonic 
bifurcation); (b) a couple of conjugate eigenvalues reaches — 1 (0 = ±7r) and leaves the 
unit circle along the real line (subharmonic bifurcation); (c) two pairs of conjugate 
eigenvalues collide at two conjugate points on the unit circle and leave it (Krein crunch 
or oscillatory bifurcation). It must be remembered that a bifurcation involving two 
eigenvalues with the same sign of the Krein signature K{d) — sign(i(^ • ^* — ' 0) 
not possible [|l[. 

The basic features of the Floquet multipliers at the anti-continuous limit for bright 
and dark breathers are the following: 

• For bright breathers, there are N — 1 Floquet multipliers corresponding to the 
oscillators at rest. They are degenerated a.t 6 = ±27ra)o/wb)- The eigenvalues at 
< 9 < TT have k > while the ones at — tt < 6 < have k < 0. li 9 — or 
TT K = 0. In addition, there are two eigenvalues at 1 -I- Oi, corresponding to the 
excited oscillator. 
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Figure 1. Floquot multipliers at zero coupling for: (loft) a bright breather and 
(right) a dark breather with the numbers of identical multipliers. 

• For dark breathers, there are two eigenvalues at 6* = ±27rwo/tfb corresponding to 
the oscillator at rest, and there are 2{N — 1) eigenvalues at 1 -|- Oi corresponding 
to the excited oscillators. 

When the coupling is switched on, the evolution of the Floquet multipliers for 
bright breathers is rather different from the case of dark breathers: 

• For a bright breather, the Floquet eigenvalues corresponding to the oscillators 
at rest lose their degeneracy and expand on two bands of eigenvalues, called 
the phonon bands. Their corresponding eigenmodes are extended. These two 
bands move on the unit circle and eventually cross each other. In this case, 
eigenvalues of different Krein signature can collide and abandon the unit circle 
through subharmonic or oscillatory bifurcations. In addition, some eigenvalues 
can abandon the bands and become localized p^ . A pair of complex conjugate 
eigenvalues can collide at 1 -I- Oi leading to a harmonic bifurcation or collide with 
the phonon band through a oscillatory bifurcation. 

• For a dark breather, the eigenvalues corresponding to the excited oscillators, with 
extended phonon eigenmodes, can either depart from the unit circle along the real 
axis (harmonic bifurcation) or move along the unit circle. In the last case, they 
can collide with the eigenvalue corresponding to the rest oscillator (with localized 
eigenmode) through a Krein crunch. Eventually, if the on-site potential is non- 
symmetric, the eigenvalues collide at — 1 + Oi, leading to a cascade of subharmonic 
bifurcations. 

The study of breather stability can be complemented by means of Aubry's band 
theory 0]. It consists in studying the linearized system (||) for ^ 0, with the 
corresponding family of Floquet operators J^e- For each operator there are 2 x N 
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Floquet multipliers. A Floquet multiplier can be written as A = exp(i6'). 9 is called the 
Floquet argument. If is real then | exp(i6')| = 1 and the corresponding eigenfunction 
of (^) is bounded and corresponds to a stability mode; if 9 is complex, it corresponds 
to an instability mode. The set of points {9, E), with 9 real, has a band structure. The 
breather is stable if there are 2 x N band intersections (including tangent points with 
their multiplicity) with the axis E = 0. The bands are reduced to the first Brillouin 
zone (— TT, tt] and are symmetric with respect to the axis 9 = 0. The fact that J-q has 
always a double +1 eigenvalue corresponding to the phase mode u{t) manifests as a 
band which is tangent to the E = axis. 

For the uncoupled system, the structure of the stable and unstable bands is 
completely explained by the theory (although it has to be calculated numerically). 
For the coupled system, such structure is expected to change in a continuous way in 
terms of the parameter e. We will use band theory, to predict the evolution of the 
eigenvalues of To, with a model parameter is changed. 

At zero coupling, the bands corresponding to the oscillators at rest can be 
analytically calculated from the equation: 

Cn + LoUn^E^,,, (11) 

they are given by E ^ ujq — uj^9'^ / {2tt)'^ . The bands corresponding to the — 1 
excited oscillators are a deformation of the band corresponding to the oscillator at 
rest (which will be called rest bands) and one of them must be tangent to E = 
axis at {9,E) = (0,0) [0. The bands are bounded from above and numbered starting 
from the top, being the 0th band the first one. If the on-site potential is soft the 
1st band will be tangent to the E = axis at (0, 0) with positive curvature. If the 
on-site potential is hard, the tangent band at (0, 0) will be the 2nd one, and will have 
negative curvature [ pT| . At the anticontinuous limit the band scheme of a bright and 
a dark breather are similar, differing only in the number of bands that corresponds 
to the oscillators at rest or to the excited oscillators {excited bands) as it is shown in 
Figures || and ^. 

The bright breather has only one band tangent to the axis E = and that band 
cannot leave this position because, when the coupling is switched on, there must be a 
tangent band corresponding to the phase mode. The situation is totally different for 
a dark breather: there are A^ — 1 bands tangent to the E = axis and A^ — 2 of them 
can move without any restriction. If some of them move upwards, the intersection 
points disappear, which implies that some Floquet arguments become complex or, 
equivalently, that some Floquet multipliers abandon the unit circle and the breather 
becomes unstable. 



5. Dark breathers with soft on site potentials 

Let us consider a model with a cubic on-site potential given by 

VM^luj^ul- ^ul (12) 

that is, (j)'{un) = — it^ in the dynamical equations (^). Figure ^ represents the Floquet 
multipliers for e 0. The left side of the figure shows the Floquet multipliers of a 
stable bright breather for a coupling e = 0.1. However, as can be seen at the right side 
of the figure, the dark breather experiences a multiple harmonic bifurcation as soon 
as the coupling is switched on. Therefore, for a lattice with a cubic on-site potential 
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Figure 2. Band scheme for bright (left) and dark (right) breathers with soft 
potential. The continuous lines correspond to the excited oscillators and are 
numbered downwards starting from zero. The dotted lines correspond to the 
oscillator at rest. 
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Figure 3. Bands scheme for a dark breather with hard potential at zero coupling. 
The continuous lines correspond to the excited oscillators and are numbered 
downwards starting from zero. The dotted lines correspond to the oscillator at 
rest. 
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Figure 4. Evolution of the Floquet multipliers with eubic on-site potential 
and attractive interaction when the coupling is switched on for: (left) the bright 
breather at e = 0.1, which is linearly stable, although reaching a possible 
bifurcation at —1; (right) the dark breather at a much smaller coupling parameter 
£ = 0.004. The breather frequency is w^, = 0.8. 



and attractive interaction, dark breathers exist but they are unstable for every value 
of the coupling. 

This problem can be investigated by means of Aubry's band analysis. It has given 
us an explanation for the previous behaviour and also a guide for modifying the model 
in order to obtain stable dark breathers. The bands at zero coupling can be seen in 
Figure ||. 

As it is explained above and shown in Figure ^, the cubic dark breather becomes 
unstable for any attractive coupling £ > through harmonic bifurcations. This is 
easily understood in terms of the band structure: N — 2 tangent bands move upwards, 
which is mathematically demonstrated in Ref. |^ . They lose the tangent points with 
the E = axis as Figure ^ shows. Therefore, in order that the breather can be stable 
all the excited bands except one have to move downwards, transforming the points 
tangent to the E = axis into intersection points. A straightforward alternative is to 
change the sign of the parameter e in (|l|). This is equivalent to using a dipole-dipole 
coupling potential W{u) = Un+i Un, i.e., the Hamiltonian can be written as: 

H = ^(^uI + ^ujIuI-^uI + ^e{un+i Un)j (13) 
with e > 0. This repulsive interaction has been used recently to model DNA-related 



dels pa, G6L 27, Eg, M, R0|. The dynamical equations with repulsive interaction 



and cubic on-site potential become: 

iln + UJ^Un - ul + e{Un-i + U„+l) = 0, (14) 



Dark breathers stability 



10 



1.5 I r 




Figure 5. Band structure at zero coupling for a cubic on-site potential. The right 
part of the figure plots the moduli of the unstable Floquet multipliers that are 
smaller than 1, being their inverses the unstable multipliers. Frequency lj^, = 0.8 
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Figure 6. Band structure of a dark breather with cubic on— site potential and 
attractive interaction at very low coupling (e = 0.004). Frequency = 0.8. 
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Figure 7. Band structure of a dark breather with cubic on— site potential and 
repulsive interaction at e = 0.015. Frequency u!^, = 0.8. 



and the linear stability equations become: 

S.n + in - 2Un £.n + e(6i-l + Cn+l) = 0. (15) 

Figure]^ displays the band structure at e = 0.015 for this system. The — 2 bands that 
are allowed to move will perform a downwards movement and therefore the breather 
is stable. 

A further increase of e leads to Krein crunches. They are caused by the mixing 
of the rest bands and the excited bands, that produces "wiggles" in the excited bands 
and gaps appears between them [ ^ . When these bands move downwards and the 
"wiggles" cross the E — axis, intersection points are lost but recovered when the 
coupling increases. It manifests as the appearance of "instability bubbles" . These 
instabilities depend on the size of the system, as will be shown in section ^. The 
system will eventually become unstable through a cascade of subharmonic bifurcations 
when the lower excited band goes below the E ^ axis, losing the intersection points 
at 6 = ±TT (see figure . 

We have obtained dark breathers for systems with cubic on-site potential and 
attractive or repulsive interaction even though dark breathers with attractive coupling 
are not stable. Figure H shows two different examples of dark breathers for this two 
types of coupling. Notice that the oscillator with small amplitude is in phase with its 
neighbours when the system is unstable (left) while is in anti-phase when it is stable 
(right). 

This type of behaviour is general for other soft on-site potentials as it can be 
obtained using, for example, the Morse potential given by 

Viun) = D{exp{-bu^)^lf. (16) 

For this system, dark breathers maintain their stability until the coupling parameter 
reaches the value e — 0.024. From this value Krein bifurcations appear. These 
bifurcation are harmless in the sense that they preserve the dark localization, but 
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Figure 8. Band structure of a dark breather with cubic on— site potential and 
repulsive interaction at e = 0.017 and lu^, = 0.8. Note the "wiggles" that appear 
when the excited bands and the rest band mix. A further increase of the coupling 
parameter e will make the excited bands fall below E = producing subharmonic 
bifurcations. 




Figure 9. Dark breathers profile for a cubic potential with attractive interaction 
(left) and with repulsive interaction(right), at £ = 0.023 . The last one is stable. 




Figure 10. Evolution of tlie Floquet eigenvalues for a Morse on— site potential and 
repulsive coupling. Left: e = 0.024, the system is still stable. Right: e = 0.035, 
the system becomes unstable due to subharmonic bifurcation and Krein crunches. 



the breather becomes quasiperiodic with a superimposed frequency to the breather 
one. For e > 0.033, the system becomes unstable due to subharmonic bifurcations. 
Figure |l^ shows this behaviour in terms of the corresponding Floquet multipliers 
for two different values of the coupling. With this potential, it has been described 
[pl that for values of the lu\, close to the linear frequency ujq (in fact, slightly above 
it), oscillatory instabilities can bring about the movement of dark breathers, when 
the corresponding eigenvalue is asymmetric and localized. It does not happen for 
frequencies far enough of luq as the ones of the order we present here. 

An interesting variation of this scheme occurs if we consider symmetric on-site 
potentials as the quartic soft potential given by 

T/K) = i^o'"n- (17) 

In this case the potential is symmetric and no subharmonic bifurcations at A = ±1 
occur because the bands are gapless at 9 = ±7r. This characteristic enlarges the 
stability range (see Figure 11). The system eventually becomes unstable (apart from 
the reentrant instabilities due to Krein crunches) through harmonic bifurcations. 



6. Dark breathers with hard on site potentials 

The band scheme of a system with hard on-site potential at e = is shown in Figure 
|. There are — 1 bands tangent to the axis E = 0. It is clear that the breather 
will remain stable provided that the tangent bands keep the intersection points when 
the degeneracy is raised with e ^ 0. In this system it means that the N — 2 bands 
that can move will perform an upwards movement. Therefore, the dark breather will 
be stable with an attractive coupling potential. We have used for the numerics the 
quartic hard on-site potential, that is 

V{u^) = ^u;lul+ ^ul (18) 

The Hamiltonian is then 
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Figure 11. Band structure for a dark breather with quartic soft on-site potential 
and repulsive coupling. Note the absence of gaps at 6 = ±tt which enlarges the 
stability range. Parameters: e = 0.018 and uj^, = 1.2. 

Figure |l^ shows that at e = 0.01 the tangent bands have moved upwards and, 
therefore, the stabihty of the system is maintained. The system loses its stabihty at 
£ = 0.022 where a harmonic bifurcation appears. The instabihty mode, shown in ^ 
with its multiphers, is an asymmetric extended one. Simulations done perturbing with 
it the dark breather give rise to an small oscillation with both sides of the chain out 
off phase, superimposed to the dark breather one, but the darkness is preserved. After 
that, there are Krein crunches due to the bands mixing with the usual properties. 

7. Finite size effects 

This final section is dedicated to comment some preliminary results relative to the 
finite size effects on instability of dark breathers. Important differences appear with 
respect to the case of bright breathers |^ . 

In the case of bright breathers, there are two kinds of size-dependent bifurcations. 
The origin of them relies on the nature of the localization of the colliding eigenvalues. 
If the colliding eigenvalues are extended, there appear "instability bubbles" , i.e. the 
Floquet eigenvalues abandon the unit circle after the collision but return afterwards. 
These bifurcations disappear when the system is infinite. Alternatively, a localized 
eigenvalues can collide with a band of extended eigenvalues. In this case, the instability 
bubbles also occur but they persist even though the system is infinite. 

In order to study the case of dark breathers, we have chosen a Morse on-site 
potential with repulsive interaction. As a result of this analysis, the collision of 
extended eigenvalues always implies a cascade of subharmonic bifurcations (see section 
^ independently of the size of the system. It is due to the fact that they are due to 
excited bands, whose Floquet arguments belong to the first Brillouin zone, losing 
intersections at 6 = tt. It is different from the case of bright breathers where this kind 




Figure 13. Zoom of the band structure for a quartic hard on— site potential with 
attractive interaction at e = 0.041 and uj\^ = 1.2. 
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Figure 14. Left: Floquet multipliers witli a quartic liard on-site potential and 
attractive interaction at £ = 0.041 and ojb = 1-2. A harmonic and a small Krein 
instabilities appear. Right: the velocities components (the position ones are zero) 
of the eigenvalue corresponding to the harmonic instability. 
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Figure 15. Zoom of the band scheme for N=21 and N=101 and a Morse potential 
(e = 0.02 and ojb = 0.8). It can be observed the appearance of wiggles when N=21, 
which are much wider when N=101 these wiggles are almost invisible because they 
occupy the width of the band. 



of instabilities are due to modes corresponding to oscillators at rest and the bands 
must be reduced to the first Brillouin zone. 

Also, an eigenvalue corresponding to an extended mode can collide with a localized 
mode. In this case, the rest (localized mode) band and the excited (extended modes) 
bands mix (see Figures |l^). There appear "wiggles" that imply the opening of gaps 
in the band scheme, which are the origin of the instability bubbles. However, when 
the size of the system is increased, the wiggles widen until they occupy almost the 
width of the quasi-continuous extended mode bands. This fact implies that, although 
the instability bubbles disappear, the Krein crunches are unavoidable (see Figure p^ . 
Nevertheless, the breathers are robust despite of the existence of these instabilities. 
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Figure 16. Moduli of the Floquet eigenvalues corresponding to the Krein 
crunches for N=21 and N=101 and a Morse potential (e = 0.02 and lu]j = 0.8). It 
can be observed the appearance of instability bubbles when N=21. However, these 
bubbles transform into an instability that persist up to a value of the coupling for 
which the breather is unstable because of the subharmonic bifurcations. 



8. Dark breathers stability in dissipative systems 

There appear some important differences in the study of the stability of dark breathers 
in dissipative systems with respect to the Hamiltonian case. 

In order to perform this study, we start from the dissipative Frenkel-Kontorova 
model 11: 

ii„ + 7u„ + V'{un) + e{2un - - = Fsmuj\,t. (20) 

where 7 is a damping parameter and F the amplitude of an external force. 

In the anticontinuous limit, there are no band tangent to the E — 0-axis 
(figure prja). It implies that, when the coupling is introduced, there are no losing 
of intersections of the bands with independence of the sign of the coupling constant 
(e). Furthermore, the band corresponding to the oscillator at rest do not mix with the 
band of the background and oscillatory instabilities do not appear. These phenomena 
are shown in figure |l^. 

9. Conclusions 

In this paper we have explored the existence and stability of dark breathers in 
one-dimensional Klein-Gordon models, for frequencies far enough from the linear 
frequencies to the DNLS approximation be justified . We have found stable dark 
breathers in several types of them. For systems with soft on-site potential, there 
are no stable dark breathers if the coupling between particles is attractive, but with 
a repulsive coupling the stability is assured for fairly high values of the coupling 
parameter e. As e increases instability bifurcations due to Krein crunches take place 
as a consequence of the band mixing between the rest oscillator band and the dark 
background. Eventually, the system experiences subharmonic bifurcations at —1 that 
make the breather unstable. If the soft on-site potential is symmetric the subharmonic 
bifurcations are avoided and the final instability is caused by harmonic bifurcations. 
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(b) 



(c) 



Figure 17. Band diagrams in tlie anticontinuous limit (top) and for a coupling 
£ = ±0.005 (bottom) for a dark breather in the Prenkel— Kontorova model. The 
parameters are: oji, = 0.2tt, 7 = 0.02 y F = 0.02. (b) corresponds to an attractive 
interaction potential, whereas (c) corresponds to a repulsive potential. The dark 
breather is stable in both cases. 



For systems with hard on-site potentials, the situation is reversed and the dark 
breathers are stable with attractive couphng and unstable with repulsive one. They 
experience harmonic bifurcations at +1 apart from the Krein crunches. An analysis of 
larger systems shows that these bifurcations persist even though the system is infinite. 
Dissipative systems are, however, stable for both types of coupling. 

We are now on the project of performing a wider study of dark breathers in 
dissipative systems. Another interesting aspect is to study the relationship of the 
vibration pattern of dark breathers with the hardness of the on-site potential and the 
type of coupling. The ansatz m„ (— l)"u„, which transforms repulsive coupling 
into attractive one and vice versa, for symmetric potentials, suggests, and it has been 
shown for the DNLS equations, and we also have checked numerically, that the stability 
conditions are reversed are low coupling, but the bifurcations for jtJb— t^ol large enough, 
highly asymmetric potentials as found in chemical and biological systems, and larger 
coupling are worth studying. 
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